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Beta Beams could address the needs of long term neutrino oscillation experiments. They can produce extremely 
pure neutrino beams through the decays of relativistic radioactive ions. The baseline scenario is described, together 
with its physics performances. Using a megaton water Cerenkov detector installed under the Frejus, Beta Beams 
could improve by a factor 200 the present limits on sin 2 2#i3 and discover leptonic CP violating effects if the CP 
phase S would be greater than 30° and 813 greater than 1°. These performances can be further improved if a 
neutrino SuperBeam generated by the SPL 4MW, 2.2 GeV, proton Linac would be fired to the same detector. 
Innovative ideas on higher and lower energy Beta Beams are also described. 



1. Introduction 

Long Baseline neutrino beams are the future 
facilities for neutrino oscillations. Only in a well 
controlled, fully optimized environment it will be 
possible to perform precision measurements of 
023 and Sm^ and to search for the still unknown 
parameters #13 , S, sign(<5m 2 ). These latter pa- 
rameters can be explored by detecting sub lead- 
ing — > v e transitions, characterized by their 
smallness given the present experimental bound, 
from the Chooz experiment pQ. 

The present generation of LBL neutrino exper- 
iments: K2K , Minos, Opera and Icarus has been 
designed to confirm the SupcrKamiokande result 
on atmospheric neutrinos through disappear- 
ance or — » v T transitions, with limited sensi- 
tivity to 013 . 

Second generation long baseline experiments, 
like the already approved T2K %, and NOz/A [3], 
will extend the sensitivity on sin 2 2^13 by more 
than one order of magnitude with respect to the 
present Chooz limit. However they will have very 
limited sensitivity to the CP phase 5 even if com- 
plemented by high sensitivity reactor experiments 

A third generation of LBL neutrino experi- 
ments will be required to start a sensitive search 
of leptonic CP violation (LCPV). These fu- 
ture experiments will push conventional neutrino 
beams to their ultimate performances (neutrino 



SuperBeams) , or will require new concepts in the 
production of neutrino beams. 

Conventional neutrino beams are generated 
by secondary particle decays, mainly pions and 
kaons, producing a multiflavour neutrino beam 
(i/ e , 77^, 77 e , besides the main neutrino compo- 
nent, ). Beam composition and fluxes are dif- 
ficult to precisely predict because of the lack of 
knowledge of secondary particle production cross 
sections (hadroproduction). 

These limitations are overcome if the neutrino 
parents can be selected, collimated and acceler- 
ated to a given energy. The neutrino beams from 
their decays would then be pure and perfectly pre- 
dictable. This can be tempted within the muon or 
a beta decaying ion lifetimes. The first approach 
brings to the Neutrino Factories [S], the second 
to the Beta Beams. 

Beta Beams have been introduced by P. Zuc- 
chclli in 2001 [Hj- The idea is to generate pure, 
well collimated and intense v e (pe) beams by pro- 
ducing, collecting, accelerating radioactive ions 
and storing them in a decay ring. The ideal candi- 
dates are chosen with a lifetime around 1 s among 
the ions that can be artificially copiously pro- 
duced. The best candidates so far are 18 Ne and 
6 He for v e and F e respectively. A baseline study 
for such a BetaBeam complex has been produced 
at CERN 

In this scenario Beta Beam neutrino energies 
are below 1 GeV and the ideal detector is a water 
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Cerenkov detector with a mass of the order of 1 
megaton, necessary to reach the needed sensitiv- 
ity. Such a detector would have excellent physics 
capabilities in its own, as fully described in ref- 
erence jS], like ultimate sensitivities for proton 
decay, supernovae neutrinos, atmospheric neutri- 
nos etc. A candidate site for this detector exists, 
located under the Frejus, at an appropriate base- 
line from CERN: 130 km. 

Liquid Argon could be an excellent alterna- 
tive to water, provided that the Icarus technology 
could be scaled to the 100 Kton scale jH]. 

2. The BetaBeam complex 

The beta-beam complex is described in jj] and 
shown schematically in figure ^ Protons are de- 




Figure 1. Schematic layout of the beta-beam 
complex. At left, the low energy part is largely 
similar to the EUPJSOL project. The central 
part (PS and SPS) uses existing facilities. At 
right, the decay ring has to be built. 



livered by a high power Linac. Beta Beam targets 
need 100 /iA proton beam, at energies between 1 
and 2 GeV. 

In case the Super Proton Linac (SPL) ^U] 
would be used, Beta Beams could be fired to the 
same detector together with a neutrino Super- 
Beam SPL is designed to deliver 2mA of 
2.2 GeV (kinetic energy) protons, in such a con- 



figuration Beta Beams would use 10% of the total 
proton intensity, leaving room to a very intense 
conventional neutrino beam. 

The targets are similar to the ones envisioned 
by EURISOL 12 : the 6 He target consists either 
of a water cooled tungsten core or of a liquid lead 
core which works as a proton to neutron converter 
surrounded by beryllium oxide |13| . aiming for 
10 15 fissions per second. 18 Ne can be produced 
by spallation reactions, in this case protons will 
directly hit a magnesium oxide target. The collec- 
tion and ionization of the ions is performed using 
the ECR technique [Hj . 

Ions are firstly accelerated to MeV/u by a Linac 
and to 300 MeV/u, in a single batch of 150 ns, 
by a rapid cycling synchrotron . 16 bunches (con- 
sisting of 2.5 10 12 ions each in the case of 6 He) 
are then accumulated into the PS, and reduced 
to 8 bunches during their acceleration to inter- 
mediate energies. The SPS will finally accelerate 
the 8 bunches to the desired energy using a new 
40 MHz RF system and the existing 200 MHz 
RF system, before ejecting them in batches of 
four 10 ns bunches into the decay ring. The SPS 
could accelerate 6 He ions at a maximum 7 value 
of 7 6 Hc = 150. 

The decay ring has a total length of 6880 m 
and straight sections of 2500 m each (36% use- 
ful length for ion decays). These dimensions are 
fixed by the need to bend 6 Heions up to 7 = 150 
using 5 T superconducting magnets. Due to the 
relativistic time dilatation, the ion lifetimes reach 
several minutes, so that stacking the ions in the 
decay ring is mandatory to get enough decays and 
hence high neutrino fluxes. The challenge is then 
to inject ions in the decay ring and merge them 
with existing high density bunches. As conven- 
tional techniques with fast elements are excluded, 
a new scheme (asymmetric merging) was specifi- 
cally conceived [TB] , 

2.1. Neutrino fluxes 

18 Ne and 6 He ions can be stored together in the 
decay ring (of course, in different bunches). Due 
to their different magnetic rigidities, these ions 
would have relativistic 7 factors in the 5 to 3 ratio, 
which is quite acceptable for the physics program. 
This will impose constraints on the lattice design 
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for the decay ring, but no impossibility has been 
identified. 

An ECR source coupled to an EURISOL target 
would produce 2 • 10 13 6 He ions per second. Tak- 
ing into account all decay losses along the acceler- 
ator complex, and estimating an overall transfer 
efficiency of 50%, one estimates that 4 ■ 10 13 ions 
would permanently reside in the final decay ring. 
That would give an antineutrino flux aimed at 
the Frejus underground laboratory of 2.1 • 10 18 
per standard year (10 7 s). 

For 18 Ne, the yield is expected to be 1.6 • 10 12 
ions in a 2 s exposing time. Due to this smaller 
yield, which could be certainly improved with 
some R&D, it was then proposed to use 3 EU- 
RISOL targets in sequence connected to the same 
ECR source. Again taking into account decay 
losses plus a 50% efficiency, this means that 4-10 13 
such ions would reside in the decay ring, giving 
rise to a neutrino flux of 0.7 • 10 18 per standard 
year. 

In the following it was supposed that the neu- 
trino flux from 18 Ne could be increased by 50% 
over the present conservative estimate, having 
room for improvements both in the cycle dura- 
tion of PS and SPS and in the 18 Ne production 
at the targets. A 40 % improvement was put on 
the 6 He generated antineutrino fluxes. 

The reference fluxes are then 2.9- 10 18 6 He use- 
ful decays/year and 1.1 ■ 10 18 Ne decays/year. 

2.2. Radiation issues 

The main losses are due to decays of He ions, 
and reach 1.2 W/m in the PS and 9 W/m in the 
decay ring. This seems manageable, although the 
use of superconducting bending magnets in the 
decay ring requires further studies. Activation is- 
sues have been recently addressed US], and show 
that the dose rate on magnets in the arcs is lim- 
ited to 2.5 mSv/h at contact after 30 days opera- 
tion and 1 day cooling. Furthermore, the induced 
radioactivity on ground water will have no impact 
on public safety. 

3. Physics reach 

Beta Beam sensitivity have been computed as- 
suming a water Cerenkov detector of 440 kt fidu- 



cial mass installed in the underground Frejus lab- 
oratory, with a 130 km baseline. Most of the re- 
sults of this section are taken from reference |17| . 

3.1. Signal and backgrounds 

The neutrino beam energy is defined by the 7 
of the parent ions in the decay ring. The energy 
optimization is a compromise between the advan- 
tages of the higher 7, as a better focusing, higher 
cross sections and higher signal efficiency and the 
advantages of the lower 7 values as the reduced 
background rates and the better match with the 
probability functions |18j . 

Given the decay ring constraint (see section 
2.1): 7( 6 He)/7( 18 Ne) = 3/5 the optimal 7 val- 
ues result to be 7( 6 He) = 60 and 7( 18 Ne) = 100. 
Fig. 121 shows the BetaBeam neutrino fluxes com- 
puted at a 130 Km baseline, together with the 
SPL Super Beam (SPL-SB) fluxes. 

The mean neutrino energies of the V e , v e beams 
are 0.24 GeV and 0.36 GeV respectively. They 
are well matched with the CERN-Frejus 130 km 
baseline. On the other hand energy resolution is 
very poor at these energies, given the influence 
of Fermi motion and other nuclear effects. Sensi- 
tivities are computed for a counting experiment 
with no energy cuts. 

A different optimization holds if each ion can be 
run separately at its optimal 7 value. In this case 
the single fluxes can be doubled just filling all the 
SPS batches with the same ion, and in 10 years 
the same integrated number of useful decays/year 
is obtained. Under this hypothesis the optimal 
choice seems to be 76fj e — 7 ls Nc = 75. 

It has been recently pointed out 19 that using 
the algorithm tools discussed in section 4, Beta 
Beam performances could be further improved by 
accelerating the ions to 7e Hc = 150 with a base- 
line of 300 km. 

The signal in a Beta Beam looking for v e — > 
oscillations would be the appearance of 
charged-current events, mainly via quasi-elastic 
interactions. These events are selected by re- 
quiring a single-ring event, the track identified 
as a muon using the standard SuperKamiokande 
identification algorithms (tightening the cut on 
the pid likelihood value), and the detection of 
the muon decay into an electron. Background 
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Figure 3. Detection efficiency for F„ (solid) and 
(dashed) events, efficiency is smaller given 
the probability for a fi~ to be absorbed before 
decaying. 
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Figure 2. Beta Beam fluxes at the Frejus location 
(130 km baseline). Also the SPL Super Beam 
and Vn fluxes are shown in the plot. 

rates and signal efficiency have been studied in a 
full simulation, using the NUANCE code [20] and 
reconstructing events in a SuperKamiokande-like 
detector [2] . The efficiency curve for and F M 
events is displayed in figure 

The Beta Beam is intrinsically free from con- 
tamination by any different neutrino flavor. How- 
ever, backgrounds can be generated by ineffi- 
ciencies in particle identification, such as mis- 
identification of pions produced in neutral cur- 
rent single-pion resonant interactions, electrons 
(positrons) mis-identified as muons, or by exter- 
nal sources such as atmospheric neutrino interac- 
tions. 

The pion background has a threshold at neu- 
trino energies of about 450 MeV, and is highly 
suppressed at the Beta Beam energies. The elec- 
tron background is almost completely suppressed 
by the request of the detection of a delayed 
Michel electron following the muon track. The 



atmospheric neutrino background can be reduced 
mainly by timing the parent ion bunches. For a 
decay ring straight sections of 2.5 km and a bunch 
length of 10 ns, which seems feasible ..T;, this back- 
ground becomes negligible 6 . Moreover, out-of- 
spill neutrino interactions can be used to normal- 
ize it to 1% accuracy level. 

Signal and background rates for a 4400 kt- 
yr exposure to 6 He and 18 Ne beams are re- 
ported in Table ^ The default values for 
the oscillation parameters are sin 2 2023 = 1, 
5ml 3 = 2.5 • 10- 3 eV 2 , sin 2 26 12 = 0.8, 5m\ 2 = 
7.1 • 10~ 5 eV 2 , sign((5™ 2 )=+l. 

3.2. Systematic errors 

Systematic errors could spoil the sensitivity of 
any experiment looking for leptonic CP viola- 
tions. To this purpose a beam without any con- 
tamination where the neutrino fluxes are known 
with great precision, like the Beta Beam, is the 
ideal facility. 

A major concern is the knowledge of neutrino 
cross-sections in an energy range where experi- 
mental data are poor and the cross sections are 
roughly proportional to -E 2 . On the other hand 
Beta Beams are the ideal place where to measure 
neutrino cross sections. Being the neutrino fluxes 
precisely known, a close detector of ~ 1 kton 
(fiducial) placed at a distance of about 1 km from 
the decay ring could directly measure the rele- 
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Table 1 

Event rates for a 4400 kt-y exposure. The signals 
are computed for 13 = 1°, 8 = 90° sign(<5m 2 ) = 
+1. "(5-oscillatcd" events indicates the difference 
between the oscillated events computed with 6 = 
90° and with 6 = 0. "Oscillated at the Chooz 
limit" events are computed for sin 2 2^13 = 0.12, 
6 = 0. 

6 He 18 Ne 
(7 = 60) (7 = 100) 



CC events (no oscillation) 
Oscillated (Chooz limit) 
Oscillated (6 = 90°, 6 13 = 3 C 
6 oscillated 
Beam background 
Detector backgrounds 



19710 
612 
44 
-9 

1 



144784 
5130 
529 
5712 


397 



vant neutrino cross sections. Furthermore the 7 
factor of the accelerated ions can be varied. In 
particular an energy scan can be initiated below 
the background production threshold, allowing a 
precise measurement of the cross sections for res- 
onant processes. 

It is estimated that a residual systematic error 
of 2% will be the final precision with which both 
the signal and the backgrounds can be evaluated. 

The 813 and 6 sensitivities are computed tak- 
ing into account a 10% error on the solar 5m 2 
and sin 2 20, and a 5% and 1% error on <5m| 3 and 
sin 2 2023 respectively, as expected from the T2K 
neutrino experiment [2]. Only the diagonal con- 
tributions of these errors are considered. 

3.3. 6*i3 1 6 sensitivities 

013 can be measured either with v e and T> e dis- 



appearance and with Vn 



appearence. The 



disappearance channels offer a cleaner extraction 
of 6*i3 from the experimental result, but are lim- 
ited by systematic errors. The comparison of the 
v e and v e disappearance measurements could also 
set limits to CPT violation effects. 

013 and 6 are so tightly coupled in the appear- 
ance channels that the sensitivity expressed for 
6 = is purely indicative. A better understand- 
ing of the sensitivity of the BetaBeam is expressed 
in the ($13 , 6) plane, having fixed all the other pa- 
rameters, as shown in Fig.0J). 
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Figure 4. 90%CL sensitivity expressed as func- 



tion of 6 for 8r 



'23 



= 2.5 • 10-W 2 . CNGS and 



J-Parc curves are taken from [231 > BNL from [21] . 
All the appearance sensitivities are computed for 
sign(5m 2 )) = +1 and 5 years of data taking. 



A search for leptonic CP violation can be per- 
formed fitting the number of muon-like events to 



the p{v e 



) and to the p(v e 



probabil- 



ities. The fit can provide the simultaneous deter- 
mination of 0i3 and 6, see figure Event rates 
are summarized in Table ^ The (0i3 ,6) parame- 
ter space where the experiment could measure a 
3cr LCPV effect is shown in figure H3 Discovery 
potential curves are drawn for 2%, 5% and 10% 
systematic errors. The 5% systematic error curve 
is roughly equivalent to the 2% curve computed 
with a detector having half the mass, showing the 
absolute need to keep systematic errors low in this 
kind of measurements. 

3.4. Parameter correlations, degeneracies 
and clones 

Correlations between 0i3 and 6 are fully ac- 
counted for, and indeed they are negligible as 
can be seen in the fits to 0i3 and 6 shown in fig- 
ure These correlations are tiny because at the 
baseline of 130 km matter effects are negligible 
and do not compete with genuine LCPV effects. 
A full computation of degeneracies, correlations 
and clones for the Beta Beam sensitivities can 
be found in [22], assuming that at the time Beta 
Beam will start no informations about the values 
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Figure 5. Fits to #13 and 5 after a 10 years run. 
Plots are shown for the ($13 ,S) values indicated 
by the stars. For the other neutrino oscillation 
parameters see the text. Lines show la, 90%, 
99% and 3a confidence levels. 



of sign(($77i 2 ) and #23 will be available. 

3.5. Synergies between the SPL- 
SuperBeam and the Beta Beam 

The Beta Beam can run with the SPL as in- 
jector, but consumes at most ~ 10% of the SPL 
protons. The fact that the average neutrino ener- 
gies of both the SuperBeam and the Beta Beam 
are below 0.5 GeV (cfr. figure [2Jl, with the Beta 
Beam tunable, offers the fascinating possibility of 
exposing the same detector to 2 x 2 beams {y^ 
and x v e and V e ) having access to CP, T and 
CPT searches in the same run. 

It is evident that the combination of the two 
beams would not only result result an increase 
of the experimental statistics, but would also of- 
fer clear advantages in the reduction of the sys- 
tematic errors and the necessary redundancy to 
firmly establish any LCP within the reach of the 
experiment. 

The CP violation sensitivities of the combined 
BetaBeam and SPL-SB experiments are shown in 
figure H3 

4. High energy scenarios 

It has been pointed out in that allowing 
for higher values of 7, that is exploiting a higher 
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Figure 6. S discovery potential at 3a (see text) 
computed for 10 years running time. Black curves 
are computed with 2% systematic errors (solid), 
5% (dashed) and 10% (dot-dashed). Dotted curve 
is the Beta Beam computed by running the two 
ions at 7 = 75 (2% systematic error). Gray solid 
curve is SPL-SuperBeam sensitivity, gray dashed 
is the Beta Beam plus SPL Superbeam sensitiv- 
ity curve, both curves are computed with a 2% 
systematic error. 



energy proton synchrotron than the SPS, perfor- 
mances of the Beta Beam could be significantly 
enhanced. In the paper three different scenar- 
ios are compared: the baseline scenario described 
in section 2 (Setup I), 76h c = 350 with a base- 
line L ~ 730 km (Setup II), and 76 Hc = 1500 
with L ~ 3000 km (Setup III). The same fluxes 
of the baseline scenario have been assumed (the 
ion currents entering the decay ring correspond 
to - 0.5 MW at 76 Hc = 350 and to - 2MW at 
7 6 Hc = 1500) and it should be noted that the 
decay ring length, linearly proportional to 76He 1 
would be of 16 (70) km for 7 6 Hc = 350(1500) if 
computed with the same parameters of the base- 
line scenario. 

The main advantages of the most promising 
Setup-II scenario would be 

• the mean neutrino energy would be E v ~ 
1.2 GeV and so water Cerenkov detectors 
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are still suitable; 

• a factor ~ 10 increase in charged current 
event rate (1.5 increase at constant accel- 
erator power) with respect to the baseline 
scenario; 

• the possibility to exploit energy spectrum 
(more powerful fits to #13 , 8) ; 

• the possibility of measuring sign(Am 2 ) 
(baseline ~ 700 km). 

On the other hand, as discussed in section 3.1, the 
background rate induced by charged pion produc- 
tion in NC events is very high at these energies. 
However the longer the pion track, the higher 
the probability it interacts in water missing the 
signature of the decay electron. A modest cut 
in the visible neutrino energy can reduce back- 
ground fractions below 3 ■ 10 -3 for an integrated 
efficiency of 30 - 50% |2S]. 

Energy reconstruction is much more demand- 
ing at 1.2 GeV, in particular the non quasi elas- 
tic fraction of CC events becomes important, but 
the methods introduced in the paper allow for an 
efficient (though indirect) true neutrino energy 
estimation from the measured charged muon mo- 
mentum. 

The 99% CL discovery potential curves in the 
three different scenarios, computed including no 
systematic errors are shown if figured 

4.1. Another high energy scenario 

In it has been introduced the case of the 
highest energy that could be accessible at LHC: 
76He = 2488. It was assumed that the same fluxes 
of the baseline scenario could be used and that 
the decay ring could be partially accommodated 
inside the LHC tunnel. The neutrino fluxes at 
the CERN-Gran Sasso baseline (732 km) would 
be so high that just counting the muons pro- 
duced by neutrino interactions in the rock, with 
a 15 x 15 m 2 counter detector, the sensitivity on 
sin 2 20i3 could be increased by about one order of 
magnitude (for 8 = 0) with respect to the baseline 
scenario. Given the neutrino energy and the base- 
line, the experiment would be strongly off-peak, 
and no information about S or sign {8m 2 ) would be 
available. 




T5 13 (dec) 

Figure 7. From reference |25| : region where 8 
can be distinguished from 8 = or 8 = 180° at 
99% CL for Setup-I (solid), Setup-II with a wa- 
ter Cerenkov detector of 400 kton (dashed) at a 
baseline of 730 km and with the same detector 
with a factor 10 smaller mass (dashed-dotted) 
and Setup-Ill (dotted) with a 40 kton tracking 
calorimeter. In all cases 10 years of running time 
are considered. 



5. Low energy scenarios 

As discussed in section 3.2, Beta Beams are 
the ideal place where to measure neutrino cross 
sections. This could allow for precision measure- 
ments in the range of interest of astrophysics (su- 
pernovae explosions, nucleosynthesis etc.): E v = 
10 — 100 MeV, corresponding to 7e Hc = 7 — 14. 
This experimental possibility is fully discussed in 

6. Conclusions 

Beta Beams are a novel concept which can pro- 
vide a very clean and powerful facility for the 
search of leptonic CP violation: a single flavour 
beam and no competition with the fake CP phe- 
nomena induced by matter effects. 

The baseline scenario has been designed to be 



8 



based on available technology with some conser- 
vative extrapolations. Its design has a very strong 
synergy with the EURISOL project aiming at 
producing high intensity radioactive beams for 
nuclear physics studies with astrophysical appli- 
cations. 

The gigantic far detector needed for these stud- 
ies has excellent physics capabilities in its own, 
i.e. to study proton decay and detect supernova 
explosions. A possible site exists near CERN at 
the right distance. 

Physics potential in neutrino oscillations would 
be a sin 2 26\3 sensitivity more than two orders 
magnitude better than the present experimental 
limit, having the distinctive feature to look for 
#13 also in v e and v e disappearance channels, and 
a 3 a discovery potential on leptonic CP violation 
for CP phase 5 values greater than 30° and for 
6*i3 greater than 1°. 

These performances can be further improved 
if a neutrino SuperBeam generated by the SPL 
4MW, 2.2 GeV, proton Linac would be fired to 
the same detector. In this configuration the two 
beams could address leptonic CP and T violation, 
and could also explore CPT violation in neutrino 
oscillations. 

Several new ideas about the baseline option op- 
timization and extensions to higher and lower en- 
ergies have been recently published, showing the 
great and partially still unexplored potential of 
the Beta Beams. 
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